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Abstract 

Ground-state of the S = 1 XXZ spin chain under the influence 
of the random magnetic fleld is studied by means of the exact- 
diagonahzation method. The S = 1/2 counterpart has been in- 
vestigated extensively so far. The easy-plane area, including the 
Haldane and the XY phases, is considered. The area suffers signif- 
icantly from the magnitude of the constituent spin. Destruction of 
the Haldane state is observed at a critical strength of the random 
fleld, which is comparable to the magnitude of the Haldane gap. 
This transition is characterized by the disappearance of the string 
order. The XY region continues until at a critical randomness, at 
which a transition of the KT-universality class occurs. These fea- 
tures are contrasted with those of the S = 1/2 counterpart. 



1 Introduction 

In quantum statistical mechanics, randomnesses still remain non- 
trivial, even though they do not introduce any competing interac- 
tions. The essence of the integer quantum hole effect, for example, 
is attributed to the randomness-driven localization-delocalization 
transition under the presence of high magnetic fleld. Although this 
description is simply of the one-body picture, it is very hard to 



deal with this problem. The effect of randomness for the quantum 
many-body system would be more complicated, and has not been 
investigated very well. In order to consider the effect of randomness 
for such the systems as the Mott insulator and the superconductor, 
for example, it is significant to take into account of the many-body 
interactions 

In one dimension, some analytic treatments are available: The 
random S = 1/2 XX Z model was investigated with the bosoniza- 
tion technique [j2j] and the real-space decimation method [Q]. The 
real-space decimation method yields comprehensive understanding 
of the random transverse-field Ising chain as well 0] . The analyti- 
cal predictions of these theories are confirmed numerically 0, §, |^ . 
In the course of the above studies, one could resort to very nice and 
special mathematical characteristics. That is, the S = 1/2 chain 
is described in the language of the bosonization formalism, while 
the transverse-field Ising chain is related to the free-fermion model 
with the Jordan- Wigner transformation. Because both the resul- 
tant models are understood very well, the effect of the randomness 
could be considered to a certain extent. 

Recently, more general systems, which are possibly missing such 
the nice characteristics as above, have been attracting considerable 
attention [|rD|, [n|, |T2|, [g, |T|, [T5], [1^, [TTl. It is suggestive that the 



real-space decimation method fails for the cases other than S = 1/2 
[T^. This is somehow consistent with the Haldane conjecture 



T^j , which states that the quantum magnet is affected by the mag- 



nitude of the spin quantitatively. According to his conjecture, the 
ground state of the S = 1 chain is massive. Hence, in this case, the 
technique []2], |T^, ^ based on the bosonization formalism fails as 
well. Few studies have been reported for the S = 1 magnet so far 
[pni , |T2], [T^, |T^, [T^ . These studies are concerned in the bond-random 
S* = 1 Heisenberg magnet. 

Here, we study the S = 1 XX Z model with the random magnetic 
field, 

L 

H = {^i^i+i + ^t^i+i + ^S-S-^i + HiS-) . (1) 

i=l 

The operators {S^} denote the S = 1 spin operators acting on the 



site i. The random field {Hi} distributes uniformly over the range 



[-V3A,x/3A]. Therefore the mean deviation is given by ^[Hf].^y — 
A. 

Under the absence of the random field A = 0, the ground state 
phase diagram of the system ([I|) is known |2T], : (It is noteworthy 
that the phase diagram is different qualitatively from that of the 
S = 1/2 chain, which is reviewed in the next section.) In the region 
A < — 1, the ground state is ferromagnetic. In the region — 1 < A < 
0, the ground state is in the XY phase. Therefore, in this area, the 
ground-state magnetic correlation decays obeying the power law, 

(2) 

The exponent 775-^1 is speculated to be described by the formula, 



TT-COS ^A . 

In the region < A < 1.2 • • •, the Haldane phase extends. In the 
Haldane phase. The ground state is massive, and the magnetic cor- 
relation decays exponentially, 

{SfS]) ~ e-l^-^'l/^. (4) 

At the isotropic point A = 1, the Haldane gap AE'tnpiet and the corre- 
lation length ^ are estimated numerically AE'tripict = 0.41049(2) 



and ^ = 6.2, respectively. The antiferromagnetic phase appears in 
the remaining region 1.2 • • • < A. Among the above four phases, 
the antiferromagnetic and the ferromagnetic phases are of rather 
classical nature. Hence, in this paper, we devote ourselves to the 
easy-plane region, namely, the Haldane and the XY phases. Both 
phases are realized precisely due to the quantum effect, and the 
influence of the randomness is unclear. 

The present paper is organized as follows. First, we review the 
results obtained for the 5=1/2 model in the next section. In Sec- 
tion 1^, our numerical results for the S = 1 chain ([I|) are presented. 
First, we focus on the transition from the Haldane phase to the lo- 
calization phase at the isotropic point A = 1. For the first time, we 
observed the onset of the localization transition from the Haldane 



phase at Ac = 0.49 ± 0.15, with the critical exponent v = 3.4 ± 2.2. 
The transition is characterized as the disappearance of the string 
correlation. In the remaining part, we concentrate on the XY re- 
gion — 1 < A < 0. We find that the XY phase persists considerably 
against the randomness, and the randomness-riven phase transition 
is of the KT universality class. In the last section, we summarize 
these results, and make a comparison between ours and those of 
S= 1/2. 



2 Review — the ground-state phase diagram 
and the criticality of the S = 1/2 XX Z chain 
with the random magnetic field 

Here, we summarize the ground-state property of the the Hamilto- 
nian ([I|) in the case of 5 = 1/2. First, we explain the phase-diagram 
without randomness [j25], |2^: the XY phase spans the whole easy- 



plane region — 1 < A < 1. In this XY phase, the critical exponent 
^5=1/2 is evaluated exactly as the function of the anisotropy A, 



TT — cos ^ A 
^5=1/2 = • (5) 

TT 

Beside this XY phase, the ferromagnetic phase appears in A < — 1, 
while the antiferromagnetic phase extends in 1 < A. 

The ground-state phase diagram in the presence of the random 
field is depicted in Fig. The diagram is determined with var- 
ious means, the bosonization technique [0], real-space decimation 
method [Q] and the exact-diagonalization method [g] . 

If a ground state is critical and the critical exponent 77 is known, 
Harris's criterion [)2^ tells whether randomness is a relevant per- 



turbation or not. The Harris criterion in the present case states 
that the random-field perturbation is relevant (irrelevant), if the 
exponent is 77 > 1/3 (77 < 1/3); and so, r/c = 1/3. Note that this 
criterion together with the formula (|5|) yields the conclusion consis- 
tent to the phase diagram shown in Fig. |T]; the XY phase is stable 
in the anisotropy range — 1 < A < —0.5, whereas it is unstable in 
the other range —0.5 < A < 1. 



The criterion, however, is concerned only with the vicinity of 
the hne A = 0. Hence, we must resort to numerical simulations 
in order to investigate the region of the finite randomness A > 0. 
The phase diagram of the full range of the randomness was explored 
by means of the numerical simulation by Runge and Zimanyi [j^. 
In our case 5 = 1, there also appears the Haldane phase, which is 
driven out of a certain criticality. Therefore, no rigorous criterion 
such as Harris's is available even in the vicinity of the pure line. 
These are the reasons why we managed the numerical simulation. 

Lastly, we mention physical interpretation of the phase diagram. 
The S = 1/2 magnet is transformed with the Jordan- Wigner trans- 
formation to the spinless-fermion system, 

n,t = E {0.5(clci+i + h.c.) + \(c\c, - l/2){cl+ic.+i - 1/2) + 0.5ff,clc.} . 

(6) 

In consequence, we see that the parameter A controls the strength 
of the repulsion among the particles, and the randomness A enters 
as the potential randomness. In general, the particles become stable 
against the randomness as the cohesive force is strengthened; the 
tendency is quite consistent with the phase diagram in Fig. 

3 Numerical results 

In this section, we investigate the ground state of the model (|I]) 
by means of the exact-diagonalization method. The results for the 
Haldane and the XY phases are presented separately in respective 
subsections. In the subsections ^TT| and each data is averaged 
over 128 and 160 samples, respectively. The diagonalization was 
performed within the sub-space of the quantum number Ef = 0. 
This sub-space is the most relevant to the ground-state: In fact, at 
A = 0, the ground state belongs to this subspace. This would not 
be violated essentially even in the presence of finite randomnesses 
A > because of the symmetry of the Hamiltonian (|I]) under the 
transformation z ^ —z. 



3.1 Random-field driven phase transition from the Hal- 
dane phase 



As is explained in Section [T], in the region < A < 1.2 • • • (A = 
0), the Haldane phase is realized. At the isotropic point A = 1, 
the magnitude of the Haldane gap becomes maximal. Here, we 
concentrate on this point A = 1. 

Some might wonder that the Haldane phase and the localiza- 
tion (random- field) phase are equivalent, because they both exhibit 
short-range magnetic correlations. Here, however, we show for the 
first time that these phases are distinguishable, in other words, a 
critical point separates these phases. 

In order to discriminate the phases, we utilized the string corre- 
lation 0, 0, 

(^stringO- -^) = [(^fe^'^n-J ^.^5;)]^^ . (7) 

Here, the bracket ( ) denotes the ground-state expectation value, 
and the symbol [ ]av stands for the random average. The string order 
develops [j3T| in the ground state of the Haldane phase, although the 
usual magnetic correlation is short-ranged; see eq. 

Because of the presence of the random field parallel to the z 
axis, the ^-component correlation is suffered from the random-field 
disturbance directly. Hence, we have evaluated the x-component 
correlation, as is defined in eq. (|7|). (In the numerical simulation, 
such a off-diagonal operator is more difficult to treat than the diag- 
onal one. This difficulty limited the maximal system size and the 
number of samples available.) 

In Fig. 0, we plotted the Binder parameter associated with 
the string correlation (|7|), 



(AC 

3 (^String/ 



TT _ -, _ \ String/ /^N 

<^ string — -L ^' 



where the string-order magnetization is given by. 



Mstring = E^fe-^-\^f. (9) 



As the system size is increased, the Binder parameter approaches 
to 2/3, if the corresponding order, namely, the string order, devel- 
ops. On the other hand, it vanishes in the disorder region, as the 
system size is enlarged. It is system-size invariant at the critical 
point. Actually, we see in Fig. g, that in the region A < 0.5 the 
Binder parameter grows through enlarging the system size, while 
it is suppressed in the opposite area 0.5 < A. In consequence, we 
see that the Haldane phase persists against the randomness up to 
Ac ~ 0.5, whereas it is disturbed by the random field in the region 
A > 0.5. 

In order to investigate the criticahty more precisely, we employed 
the finite-size-scaling technique. The finite-size-scaling theory tells 
that a dimensionless quantity such as the Binder parameter should 
be described by the formula, 

U{L) = U{L/0. (10) 

The correlation length diverges in the vicinity of the critical point 
Ac in the form, 

^ = |A-Acr^ (11) 

where u is the correlation-length critical exponent. As a conse- 
quence, the [(A - Ac)Li/^]-t/ plots should fo rm an universal curve 
irrespective to the system size L. In other words, the scaling param- 
eters, such as Ac and v, are chosen so that the plots, the so-called 
scaling plots, align all along a curve. The degree of the alignment 
is measured in the way explained in Appendix |^. 

The scaling plot is presented in Fig. 0. The plot determined 
the critical randomness Ac = 0.49 ± 0.15 and the critical exponent 
u = 3.4 ± 2.2. We find that the magnitude of the critical point is 
comparable with the magnetic excitation gap AE'tripict = 0.41049(2) 
21] at the condition A = and A = 1. 



3.2 Random-field-driven transition from the XY phase 



Here we investigate the XY phase — 1 < A < 0. The XY phase is 
characterized by the presence of the spin stiffness 



Pi 



= L 



- 



(12) 



The symbol Eo{9) denotes the ground state energy, and the angle 
6 is associated with the gauge twist of the boundary condition; 
Si Si e^^S^S^ . In other words, the spin stiffness measures the 
long-range coherence of the gauge symmetry, and thus it is related 
to the superfluid density in the context of the electron system. The 
spin stiffness should vanish in the localization (random-field) phase, 
because the global coherence is destructed there. 

In Fig. g, we plotted the spin stiffness against the randomness. 
The stiffness is suppressed, as the randomness is strengthened. In 
the weak-random region A < 1, the stiffness hardly changes with 
respect to the system size. Hence, we observe in this area, the 
stiffness remains finite in the thermodynamic limit; that is, in this 
region, the XY phase is realized. In the area A > 1, however, the 
stiffness vanishes rapidly through enlarging the system size. Each 
spin is enforced to point towards the random-field direction; the 
localization phase is realized there. To summarize, we see that the 
XY phase persists up to the critical randomness Ac ~ 1, and beyond 
the threshold the ground state becomes localized. 

We explore the criticality with the finite-size-scaling method. Be- 
cause the stiffness is system-size invariant at the critical point, the 
stiffness is described by the formula, 

P. = Ps{L/0. (13) 

(Note that it is identical with the form for the Binder parameter 
(pup.) The correlation length ^ would obey the KT type form, 

because a finite critical region extends as is shown above. 

The scaling data for A = —0.5 is presented in Fig. §. This plot 
yields the estimates Ac = 0.85 ±0.30 and A ~ 0.65. (For the details 
of the analysis, refer to Appendix 0.) 

In order to confirm the above result, we show another analysis for 
determining the transition point. As is reviewed in Section |2|, the 
Harris criterion insists that the correlation decays with the exponent 
r/c = 1/3 at the localization-delocalization transition point. In order 
to adopt this criterion, we plotted the log-log plot of the square of 



the normalized staggered magnetization, 

MijL' = l/L'[j:{-iyS^ ^1/L\ (15) 

for various randomnesses in Fig. p. We observe that the above 
criterion yields the estimate Ac ~ 1. This conclusion is consistent 
with the former estimate Ac = 0.85 ± 0.30. The agreement would 
be satisfactory. 

Secondly, we present the results at the condition A = 0; see Figs. 
|7|-0. The results resemble those of A = —0.5: The scahng analysis 
for the spin stiffness yields the critical point Ac = 0.87 ± 0.30 (see 
Fig. g), whereas the Harris criterion suggests that Ac ~ 1 (see Fig. 
|9|). Again, both estimates are in good agreement. 

Lastly, we show the results for A = —1; see Figs. [ItlHT2] . (We 



found that in the range A > 0.6 of L = 14, diagonalization con- 
vergence is hardly achieved for some random samples. Hence, the 
corresponding data are missing in these figures.) We observe in Fig. 
[TD| that solely at A = 0, XY state is realized, and for the region 
other than that, the localization phase is realized. In fact, we find 
that the data is scaled well under the assumption of the formula ( P~r|) ; 
see Fig. |TT|. The analysis gives the criticality Ac = —0.13 ±0.10 and 
u = 0.70 ±0.10. The negative value of the randomness is apparently 
non-sense. Considering a certain statistical error, we suspect that 
the transition point is located at A ~ 0. The Harris criterion, see 
Fig. [T2], confirms this conclusion: In the figure, it is shown clearly 
that the correlation is suppressed owing to the randomness signifi- 
cantly. That is, all the data are curved convexly. This fact suggests 
that the correlation decays exponentially even for the infinitesimal 
randomness. 



4 Summary and discussions 

We have investigated the S = 1 XX Z chain with the random mag- 
netic field (|I]) by means of the exact-diagonalization method. To 
summarize our numerical results presented in Section 0, we depicted 
a drawing of the ground-state phase diagram as in Fig. |T3|. This 



diagram should be contrasted to that of the 5 = 1/2 model shown 
in Fig. |I]. 

First, we discuss the Haldane phase (A = 1). Although both the 
Haldane and the random-field phase are magnetically disordered, 
these phases are separated distinctively by the phase boundary. The 
Haldane phase is characterized by the presence of the string order 
(|7|), whereas the localization phase is missing the order. We found 
that the phase boundary is located at the random-field strength 
Ac = 0.49 ±0.15. It is noteworthy that this critical value is compa- 
rable with the magnitude of the Haldane gap, AE'tnpict = 0.41049(2) 
2^] . Moreover, we obtained the critical exponent u = 3.4 ±2.2 with 



the finite-size-scaling analysis. At present, little is reported on the 
ground-state criticality driven by randomness on related systems: 
In the case oi S = 1/2 at the isotropic point A = 1, the localization- 
delocalization transition occurs at Ac = with the exponent v = 1 
[§, @]. Our universality class would not be identical to this, and 
possibly unique. 

Lastly, we mention the XY phase (—1 < A < 0). The phase 
boundary is determined both by the scaling plot of the spin stiff- 
ness and the Harris criterion on the correlation-decay rate. Both 
estimates are in good agreement. We observed that the whole re- 
gion — 1 < A < is stable against the randomness, and the phase 
boundary extends up to the random strength A ~ 0.9. Note that 
in the case oi S = 1/2, the XY phase is stable only in the range 
— 1 < A < —0.5, and furthermore the maximal critical randomness is 
~ 0.12 at most. These maximal critical randomnesses are far apart. 
Even if we normalize the energy unit in terms of the magnitude of 
the respective spins, the discrepancy between these estimates is not 
resolved completely. We conjecture that through enlarging the mag- 
nitude of the constituent spin, classical nature emerges very rapidly. 
That is to say, for a large-spin chain, the ground-state magnetism 
is very stable in respect to the randomness. This consequence is 
quite reasonable, because classical magnetic order survives even in 
the presence of the random field at the ground state. 

We have not dealt with the region a bit exceeding the point A = 0. 
At this point (A = 0), the KT transition takes place; The energy 



gap starts to open very slowly. Hence, severe correction to the 
finite-size scaling emerges. This difficulty is remained to be solved 
for the future. 
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A Details of the present scaling analyses 

We explain the details of our finite-size-scaling analyses, which we 
managed in Section ^ in order to estimate the transition point Ac, 
the exponent u in eq. (|TT]), and the coefficient A in eq. ( pT^) . We 
adjusted these scafing parameters so that the scaled data, such as 
those shown in Figs. and|5], form a curve irrespective of the system 
sizes. In order to see quantitatively to what extent these data align, 
we employ the "local linearity function" S defined by Kawashima 
and Ito Suppose a set of the data points {{xi^yi)} with the 

errorbar {di{= 5yi)}^ which we number so that xi < Xi+i may hold 
for i = 1, 2, • • • , n — 1. For this data set, the local- linearity function 
is defined as 

n-l 

S = w{x^, yi, d,\x,^i, y^^i, d^^i, x.+i, y^+i, d^+i). (16) 

i=2 

The quantity w{xj, yj, dj\x^, y^.di.Xk.yk, 4) is given by 

» = f^y. (17) 



where 



and 



A 



y = i-Loj 



Xk Xi 
2 

-Xk Xj^ / \Xk Xj 



^' = 4+(!^-^^0 • (19) 



In other words, the numerator yj — y denotes the deviation of the 
point {xj^yj) from the line passing two points {xi^yi) and {xk^yk), 
and the denominator A stands for the statistical error of {yi — y). 
And so, w = {{yi — y)/A)'^ shows a degree to what extent these three 
points align. The advantage in the above analysis is as follows: 
In the conventional least square fitting, we need to assume some 
particular fitting function. The assumption which function we use 
causes systematic error. Note that in the present analysis, we do 
not have to assume any fitting functions. 

The scaling parameters are determined so as to minimize the 
local-linearity function. The error margin of the scaling parameter 
is hard to estimate. It is concerned with both the statistical error 
and the correction to the finite size scaling. The former error can be 
estimated through considering the statistical error of the function S. 
This function has a relative error of the order 1 / — 2. (The local- 
linearity function is of the order (n — 2), and the error is given by 
y/n — 2.) The correction to the finite-size scaling, on the contrary, 
is quite difficult to determine. In the present study, in particular, 
the KT transition occurs. In that case, in general, the correction to 
the finite-size scaling emerges severly. Here, we consider only the 
statistical error in order to estimate the error margins of the scaling 
parameters. As the number of the data points n is increased, the 
statistical error of S would be reduced. The corrections to the finite- 
size scaling might increase instead. In the present analyses, we used 
twenty data in the vicinity of the transition point Ac. An example 
of the plot S is shown in Fig. [T^ We observe the minimum at 
ly ~ 1.5. Taking into account of the statistical error, we estimated 
the critical exponent as u = 3.4 ± 2.2. 
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Figure 1: The ground-state phase diagram of the S = 1/2 XX Z model (|1]) 
under the presence of the random magnetic field A. Here, we have referred to 
the figure appearing in the article 0. 




Figure 2: Plot of the Binder parameter (^) associated with the string order (|^ 
with the randomness A varied at the isotropic point A = 1. The system-size- 
invariant point indicates the location of the critical point. 
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Figure 3: The scaling plot of the data shown in Fig. |[ This plot yields the 
estimates Ae = 0.49 ± 0.15 and u = 3A± 2.2. 




Figure 4: Plot of the spin stiffness (0) against the randomness at A = 



-0.5. 
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Figure 5: Scaling plot of the data shown in Fig. ^. This plot yields the estimate 
Ac = 0.85 ±0.30. 




Figure 6: The square of the staggered magnetic moment (^) is plotted for 
various randomnesses at A = —0.5. The dashed slope shows the decay rate 
r]c = 1/3, which is proposed for the criterion for the localization-delocalization 
transition. 




Figure 7: Plot of the spin stiffness (|T2D against the randomness at A = 0. 
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Figure 8: Scaling plot of the data shown in Fig. ^. This plot yields the estimate 
Ac = 0.87 ±0.30. 




Figure 9: The square of the staggered magnetic moment (^) is plotted for 
various randomnesses at A = 0. The dashed slope show the decay rate r/c = 1/3. 
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Figure 11: Scaling plot of the data shown in Fig. 0. This plot yields the 
estimate Ac = -0.13 ± 0.10. 




Figure 12: The square of the staggered magnetic moment (|I5]) is plotted for 
various randomnesses at A = — 1. The dashed slope show the decay rate r/c = 1/3. 
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Figure 13: Phase diagram of the S = 1 XX Z model with the random magnetic 
field; the Hamiltonian is given by eq. (|l|). 




Figure 14: Local linearity function S (|I6|) for the scaling data shown in Fig. |^ 
is plotted for various trial values of u with Ac = 0.49 fixed. The location of the 
minimum yields the estimate of the exponent ~ 1.5. 



